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INTRODUCTION 


This research is part of the recent AFOSR initiative in Biomimetics, which has the broad 
overall goal of seeking to elucidate and understand the salient characteristics of naturally occuring 
biological composite materials with an aim toward mimicking these features in devising improved 
man-made composite materials. The biological composite of interest in this effort is bone tissue, 
or more specifically Haversian compact bone tissue. The main objective of the project is to 
develop a deterministic understanding of the relationships between microstructural details and 
macroscopic mechanical behavior. The unique fiber/matrix interphase material, which is 
distinctive of this type of bone tissue, is of particular interest because of the recognized 
significance of interfaces in developing new higher performance composites. The methodolog>- 
combines composite micromechanics modeling techniques to analytically model structure/property 
relationships with companion experimental testing of tissue samples to evaluate and validate the 
modeling. 

The specific aims of the first year of the project were: 

■ To develop and evaluate a microstructural model of Haversian cortical bone tissue 

which includes major microstructural features (e g., Haversian porosity, a distinct 
cement line interphase material, interface conditions between the cement line and the 
osteon and interstitial bone matrix). 

■ To establish the appropriate experimental protocol and procedures for measuring 

macroscopic mechanical properties. (Preliminary testing will include both direct 
mechanical tensile testing of bone specimens as well as ultrasound analysis.} 

■ To investigate various options for quantitative characterization of microstructural details 

such as porosity (Haversian, lacunal and other sources), osteon density, cement line 
density, percent mineralization, and overall wet and dry densities. 


RESEARCH DESCRIPTION 


Backoround 

The particular type of bone tissue being studied in this research is Haversian compact 
bone tissue because its microstructure resembles that of a fiber reinforced composite (Fig. 1). 
The fibers in this case, however, are hollow, layered, and surrounded by a distinct "interphase" 
fiber/matrix material, known as the cement line. The hollow fibers are called Haversian systems, 
or secondary osteons, and or not perfectly straight but in fact spiral slightly and branch. The 
"matrix" material is interstitial bone lying between complete osteons. Interstitial bone is actually 
composed of remnants of resorbed osteons and is therefore similar in composition and properties 
to the "fiber" component. The cement line is an amorphous substance that is typically considered 
to be very compliant or even viscous. Compact bone tissue is the more dense material found in 
the mid-shaft (diaphysis) of major load-carrying bones of humans and many mammals, as 
opposed to the less dense 'spongy' trabecular bone in the bulbous ends of these bones. 
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Micromechanics Modeling 

Finite element unit ceil analysis . A two-dimensional finite element model assuming 
hexagonal packing was developed and results compared with a previous study assuming square 
packing (Hogan, 1992). The idealized microstructure and unit cell are presented in Fig. 2. A 
typical mesh generated for this case is shown in Fig. 3. The solid modeling package PATRAN 
was used in conjunction with the ABAQUS finite element solver to analyze this model. The 
mechanical properties of each constituent are summarized in Table 1. Table 2 compares the 
hexagonal array results with those generated using a square array. Note that Ej, is the 
longitudinal, or axial, modulus whereas E 22 and E 33 are moduli in the transverse directions. The 
Haversian porosity in this case is 3%. (The results are essentially the same for hexagonal and 
square packing with the exception of a slight discrepancy for £22 ) In addition, hexagonal 
packing offers the advantage of being able to model higher fiber volume percentages, so it has 
been chosen as the method for all subsequent finite element unit cell modeling. 


unit cell 


mms 



osteon 


cement line 


Haversian canal 


Figure 2. Idealized microstructure of 
Haversian compact bone. 


Figure 3. Finite element mesh analyzed. 







Table 1 - Constituent properties used in micromechanics model 


Constituent 

E (GPa) 

V 

Osteon 

12.0 

0.3 

Cement line 

6.0 

0.4 

Interstitial bone 

15.0 

0.3 


Table 2 - Macroscopic properties for hexagonal vs. square packing geometry 


Array 

Eu 

E22 

E 33 

V 12 

V 23 

V 21 

Square 

12.3 

11.7 

11.7 

0.30 

0.30 

0.29 

Hexagonal 

12.3 

11.3 

11.3 

0.30 

0.31 

0.28 


Note: elastic moduli are in GPa 


The Haversian porosity is a parameter which greatly influences the macroscopic material 
properties of bone and is considered an important variable in micromechanics modeling. In the 
finite element model, the Haversian porosity is altered by changing the diameter of the Haversian 
canal, i.e. the inner radius of the hollow fiber. Porosity variations from 3% to 9%, which are 
reasonable values found in the literature, have been studied. Table 3 details the influence of 
porosity on transverse and longitudinal moduli as well as Poisson's ratio. The results are as 
expected: the moduli generally decrease with increasing porosity, and the transverse moduli are 
more sensitive to porosity variations. 


Table 3 - Effect of porosity on macroscopic properties (Square array) 


Porosity (%) 

E„ 

E 22 

V 12 

'^23 

V 21 

3 

12.3 

11.7 

0.30 

0.30 

0.29 

5 

12.2 

11.1 

0.30 

0.29 

0.28 

7 

11.9 

10.6 

0.30 

0.29 

0.27 

9 

11.6 

10.1 

0.30 

0.28 

0.26 


Note, elastic moduli are in GPa 
























































Separating the osteon from the interstitial bone is a somewhat mysterious substance known 
as the cement line. The exact mechanical properties (modulus and Poisson's ratio) of the cement 
line are not known, even in general, from experimental testing and can only be inferred from 
indirect testing and from the observed composition of the substance. It is generally believed that 
the cement line is more compliant (Burr et al., 1988) than either the osteon or the interstitial bone 
and that it may in fact be viscous in nature (Lakes and Saha, 1979). Moreover, the precise role 
of the cement line in determining the macroscopic mechanical behavior of compact is not known 
and much debated. Accordingly, cement line modulus values ranging widely (6.0 to 0.12 GPa, 
or, 50% to 1% of the osteon modulus) have been modeled and the resultant macroscopic 
properties documented. The thickness of the cement line was also varied from 1.4nm to S.Ofim 
since it is reported to be between l|i.m and 5(j.m (Martin and Burr, 1989). The effects of these 
changes on macroscopic properties are shown in Table 4 (1.4p.m cement line) and Table 5. As 
expected, the cement line affects the transverse moduli much more than the longitudinal moduli. 
The transverse moduli are particularly sensitive to cement line thickness. Another interesting 
observation is that these results suggest that the cement line modulus is most likely not as low 
as 0.12 GPa (i.e. 1% of the osteon modulus), because the macroscopic transverse modulus is 
lower (6.38 & 3.66 GPa) than generally observed experimentally. 


Table 4 - Effect of cement line modulus on macroscopic properties 



E„ (GPa) 

El,(GPa) 

E22(GPa) 

V23 

V 21 

6.0 (50%)* 

12.3 

11.5 

0.31 

0.28 

3.6 (30%) 

12.1 

11.3 

0.31 

0.28 

1.2 (10%) 

11.9 

10.4 

0.33 

0.26 

0.12 (1%) 

11.6 

6.38 

0.36 

0.16 



* percent of osteon modulus 


Table 5 - Effect of cement line thickness on macroscopic properties 



Ec. 

Eu ( 

GPa) 

E 22 ( 

GPa) 

V. 

23 

1 

V 

2i 

(GPa) ■ 

1.4pin 

5.0)Xm 

1.4pm 

5.0)im 

1.4pm 

5.0pm 

1.4pm 

5.0pm 

6.0 (50%)* 

12.3 

11.8 

11.5 

10.9 

0.31 

0.32 

0.28 

0.29 

3.6 (30%) 

12.1 

11.6 

11.3 

10.3 

0.31 

0.33 

0.28 

0.27 

1.2 (10%) 

11.9 

11.4 

10.4 

8.53 

0.33 

0.36 

0.26 

0.23 

0.12 (1 %) 

11.6 

11.3 

6.38 

3.66 

0.36 

0.35 

0.16 

0.10 
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Mori-Tanaka effective medium approach . The Mori-Tanaka method of estimating the 
overall mechanical properties of a fiber reinforced composite has been extended to allow for 
porosity by Zhao, et al. (1989). Even though experimental data is limited, reasonable predictions 
were found for variations in mechanical properties as a function of porosity. The simplest case 
to consider as a model for Haversian compact bone is that of 'needle-shaped' voids representing 
the Haversian canal. That is, the voids are long and slender, parallel to one another, and aligned 
with the overall longitudinal axis. The variation in longitudinal modulus with porosity is linear 
in this case, however, whereas experimental data suggests a power law relation (Schaffler and 
Burr, 1988). 

More plausible comparisons with experimental data are nevertheless possible by using the 
Mori-Tanaka method for 3-D ramdomly oriented spheroidal voids with appropriate aspect ratios. 
Results are shown in Figure 4 for several aspect ratios. The ratio of the longitudinal modulus 
to the modulus of the solid phase (E,,/Eo) is plotted versus porosity. The aspect ratio (a) of 10 
represents a prolate spheroidal void and the other three aspect ratios (.10, .05, .02) represent 
various oblate spheroidal voids. Prolate voids are more 'rod-like', while the oblate voids are more 
'disk-like'. The power-law expression used to fit the experimental data is given by: 

£ = 33.88(1 {inGPa) 

The curves in Fig. 4 that best fit the data are for oblate voids. This certainly does not describe 
the shape of the Haversian canal, but the experimental data is for bovine bone, which contains 
both Haversian and plexiform types of compact bone. No experimental data is currently available 
for modulus versus porosity for exclusively Haversian compact. The experiments to be conducted 
during the second year of this research project wnll provide this data. 



Figure 4. Longitudinal modulus vs. porosity. 
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A more realistic approach would be to model the solid portion of the bone microstructure 
as distinct components (osteon, cement line, & interstitial bone) rather than as a single phase. 
The work of Dasgupta and Bhandarkar 0 992) considers the Mori-Tanaka technique for the case 
of a fiber surrounded by one or more 'interphases' embedded in a matrix material. This model 
has been exploited for the case of Haversian cortical bone by regarding the Haversian canal as 
the fiber with the osteon and cement line as two interphase materials. Results from this approach 
are compared with those of the finite element unit cell approach in Fig. 5. The ratios of the 
longitudinal to transverse moduli are plotted vs. cement line modulus for two cement line 
thicknesses. The methods two match well for higher cement line moduli, but the Mori-Tanaka 
approach seriously underestimates the transverse modulus otherwise. 


Experimental Testino 

Specimens and preliminary testing . A working arrangement has been developed between 
the project and several laboratories in order to establish an experimental protocol for specimen 
acquistion, storage, machining, and testing. Fresh equine bone samples are obtained from the 
necropsy lab in the Department of Veterinary Physiology and Pharmacology at Texas A&M 
University. Sections of the middle portion of the cannon bone are cleaned and cut to lengths of 
approximately 2'/j inches. The specimens are kept frozen until further cutting. They are then 
rough cut into quarters and placed in Ringer’s solution and refrigerated. An Isomet diamond 
blade wafering saw is used next to cut rectangular slabs (approx. 2'x%"xW). These slabs are then 
submitted to the Testing, Maching, and Repair Facility of the Texas Engineering Experiment 
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Figure 5. Comparison of micromechanics techniques. 


9 




















Station for machining into standard "dogbone" shaped tensile specimens (see Fig. 6). The 
machining is done at low speed and constant irrigation with distilled water. Samples of human 
bone also been acquired from the University of Texas Medical Branch at Galveston under the 
subcontract. Two sections of tibia and femur, and one radius and ulna have been received. 

Once machined, each specimen is stored in a separate vial of Ringer's solution and 
refrigerated until tested. Testing must follow in less than 90 days or the specimens should be 
refrozen. Custom extensometers have been designed and fabricated for measuring the extension 
of the gage section. During preliminary testing, strain gages are also being placed on the front 
and back of the specimens. The gages are 'T gages’ to allow determination of both axial and 
transverse strains. This permits calculation of longitudinal and transverse elastic moduli as well 
as Poisson's ratio. Several Instron and MTS testing machines are available for testing the 
specimens. Four specimens (2 equine and 2 human) have been prepared thus far, and the 2 
equine speciintris have been tested. Preliminary testing is continuing in order to further 
streamline the protocol and evaluate the cost/benefit of including strain gages on every specimen. 



Figure 6. Specimen dimensions (all in inches). 
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Image analysis of microstructure . The compositional and microstructural characteristics 
will be determined quantitatively for each specimen tested. Following testing, a short piece of 
the gage section will be cut near the fracture surface and a thin section removed for 
microstructural analysis. The remainder will be used to determine the wet and dry density, and 
the percent mineralization (ash content). Facilities for computerized image analysis of the 
specimen microstructure have been obtained by the subcontractor. A Optimas image analysis 
system was purchased and procedures have been developed for automated, or semi-automated, 
determination of percent porosity, percent Haversian area, cement line density, and predominant 
collagen fiber orientation within osteons. This data will be used to correlate with mechanical 
properties and thereby evaluate the efficacy of the micromechanics modeling efforts. 

Work has also been initiated during the first year to characterize additional features of the 
porosity. The main source of porosity is considered to be the void space of the Haversian canals 
lying at the center of each osteon. These canals are generally oriented parallel with the main 
longitudinal axis of the bone, but considerable misalignment and branching actually exists. Thus, 
a more complete characterization of porosity would include some measure of the architecture of 
the canals in addition to the mere volume percentage. To this end, 3-D image reconstruction of 
Haversian systems is being pursued. This procedure involves sectioning of a portion of the 
specimen near the fracture surface and tracing the paths of several osteons along the length of 
the specimen. The 3-D computer image can be used to calculate quantitative parameters such 
as misalignment angle, degree of taper of Haversian canal, and number of branches. 
Collaboration has been initiated with researchers in the Department of Veterinary Anatomy and 
Public Health at Texas A&M University to develop and implement the image reconstruction 
procedures. Microscope slides containing a series of sections spaced 5pm apart have been made 
from specimens used in the preliminary testing. Images from the slides have been successfully 
captured to computer and imported into image analysis software packages. The same procedure 
will be repeated with human bone specimens before actual reconstruction will be attempted. 


Dissemination of Results 


The following presentations and publications have resulted from work done during the 

first year of the project; 

2-D unit cell modeling - De Frese, R.J. and Hogan, H.A., "The Influence of Packing Geometry 
and Fiber/Matrix Interface Conditions on Composite Micromechanics Modeling of 
Compact Bone Tissue", 1992 ASME Applied Mechanics, Materials and Aerospace 
Summer Meeting, April 28-May I, 1992. Scottsdale, AZ. 

Porosity modeling - Hogan, H.A. and De Frese, R.J., "A Comparison of Methods for Including 
Porosity in Micromechanics Modeling of Compact Bone Tissue", 1992 ASME Winter 
Annual Meeting (AERO 10. Biomimetics), November 8-13, 1992, Anaheim, CA ^ 

Cement line influence - De Frese, R.J. and Hogan, H.A., "A Micromechanics Study of the 
Influence of the Cement Line on the Anisotropy of Compact Bone Tissue", to be 
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presented at the 1993 ASME/AICHE/ASCE Summer Bioengineering Conference, June 25- 
29, 1993, Breckenridge, CO.^ 

Micromechanical modeling of porosity - Hogan, H.A. and De Frese, R.J., "Theoretical Modeling 
of the Effects of Porosity on the Mechanical Properties of Compact Bone Tissue", to be 
presented at the 1993 International Society of Biomechanics Congress, July 4-8, 1993, 
Paris. France.^ 

^ copies of manuscripts are included in the Appendix 


FUTURE WORK 


The goals for the second year of the project can be summarized as follows; 

1. To continue and complete mechanical testing of human, equine, and possibly bovine, bone 

tissue for the purpose of establishing a statistically meaningful data base of material 
properties which can be used to verify micromechanical models. The target number of 
specimens is 10 to 20 of each specimen type. 

2. To quantify the compositional and microstructural characteristics of each specimen tested and 

correlate this data with macroscopic mechanical properties. 

3. To improve the finite element micromechanics modeling by incorporating interface slip 

between the cement line and the surrounding bone (osteon or interstitial matrix), modeling 
the cement line as a viscoelastic interphase material, and radially varying the material 
properties of the osteon (treating it as a composite layup). 

4. To continue to investigate other micromechanics approaches such as Mori-Tanaka, self- 

consistent schemes, and homogenization theory, the latter of which has been used by 
others (Crolet et al., 1993) recently for heirarchical modeling of cortical bone 
micromechanics. 

5. To initiate ultrasonic testing of the bone specimens with the goal of obtaining transverse 

elastic moduli that can then also be correlated with microstructural characteristics. 

6. To use the micromechanics modeling to study the mechanical role of the cement line and 

assess its significance and potential as a major biomimicking feature of Haversian cortical 
bone tissue. 
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Department of Mechanical Engineering 
Texas A&M University 
College Station, Texas 


INTRODUCTION 

The next generation of advanced composite materials will have to meet increasingly 
severe and demanding performance requirements. Their successful development will 
therefore require even more rigorous and innovative development procedures. 'Tailoring' 
microstructural details, such as constituent composition and arrangement, in order to achieve 
a particular set of macroscopic mechanical properties is one such promising approach. The 
ultimate success of these and similar efforts, however, depends largely on an in-depth and 
deterministic knowledge of structure/property relations and their cause and effect 
mechanisms. In pursuing a greater understanding of these basic principles, a vast resource 
of valuable information is available through the world of naturally occurring biological 
composite materials. These highly complex, extremely efficient, and functionally adaptive 
‘natural’ composites include such materials as the intricate layered shells of many marine 
animals, the fibrous structure of wood and plants, and the highly mineralized vertebrate 
skeletal system. Tlie focus of the present work is on a particular type of bone tissue. 
Haversian compact bone, which comprises the dense bone in the mid-shaft region of the 
long bones of the legs and arms. More specifically, the emphasis is on developing composite 
material micromechanics modeling and analysis techniques as a mechanism for gaining 
insight into the fundamental relationships between bone microstructure and macroscopic 
mechanical behavior. 

A section of compact bone in the tubular mid-shaft portion of a femur is shown in 
Fig. l(a.). The major microstructural features are included, although somewhat idealized. 
This bone tissue can be considered to resemble a fiber reinforced composite in many 
respects but with some important exceptions. In adult human bone, secondary osteons (also 
called Haversian systems) predominate and are considered to represent the ‘fiber’ 
component. The fibers in this case however are hollow, less regularly spaced, and actually 
composed of several concentric layers of bone material (called lamellar bone). Osteons 
tend to be roughly aligned with the main axis of the long bone but also spiral somewhat at 
an angle with this axis. Haversian canals form the hollow centers of osteons, whereas 
Volkmann’s canals tend to traverse radially (or transversely) through the osteons connecting 
the Haversian canals of adjacent osteons. The ‘matrix’ is also made up of lamellar bone and 
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Several features of Haversian compact bone make it a particularly attractive 
candidate material for ‘biomimicking’, i.e. trying to identify fundamental charaaeristics of 
‘biological composites’ that can be exploited in developing new man-made composite 
materials. The mechanical properties of the fiber and matrix components in compact bone 
arc of the same order of magnitude since the materials are quite similar (osteons and 
interstitial bone). A similar situation exists for many ceramic and metal-matrix composites. 
Furthermore, a common toughening mechanism for such composites is the creation of weak 
fiber/matrix interface, whereas this already exists to a certain extent in compact bone. The 
mechanical significance of the cement line material of compact bone is not fully understood 
but is believed to play a similar role [3]. Many consider the cement line to be much more 
compliant than the surrounding bone (4,5] and perhaps even viscous-like [ 6 ]. Thus, it acts 
as a crack arrestor and probably plays a greater role in determining the dynamic behavior 
of compact bone. The osteons in compact bone are anisotropic due to the orientation of 
collagen fibers in the concentric layers of lamellar bone, which can be considered to be a 
functionally graded, or tailored, fiber component. Composite micromechanics modeling is 
a promising approach for quantifying these effects, but a confounding factor that must be 
included is porosity, or void space. The mechanical effects of porosity must be adequately 
accounted for before the advantages of micromechanics modeling can be fully realized for 
studying further structure/property relations. The potential medical benefits of a valid 
micromechanics model for compact bone should not be overlooked as well. Such a model 
would permit prediction of changes in mechanical properties due to mechanical adaptation 
or pathological conditions such as osteoporosis. 


ANALYSIS 

Einile. Element UniLCell Approach 

A composite material micromechanics model for compact bone has recently been 
developed and studied (7]. A quarter-fiber finite element micromechanics model for fiber 
reinforced composites was adapted and modified for this purpose. Alternate configurations 
have also been investigated [ 8 ]. Predictions of macroscopic mechanical properties are 
reasonable but somewhat difficult to judge since the mechanical properties of the 
constituents are not well established. The components of the quarter-fiber unit cell are 
shown in Fig. 2. Porosity due to the Haversian canal is represented explicitly by modeling 
the fiber component as hollow. The interstitial bone and cement line were treated as 
isotropic materials with elastic moduli of I5GPa and 6 GPa, respectively. The osteon was 
assigned a longitudinal modulus (f,,) of 24GPa and a transverse modulus (£, 2 =£ 33 ) of 
12GPa. Poisson’s ratios were taken to be 0.4 for the cement line and 0.3 for the osteon and 
interstitial bone components. 

The unit cell depicted in Fig. 2 was modeled using 2-D continuum finite elements. 
A generalized plane strain formulation was employed to permit out-of-plane loading and 
deformation. The displacements of nodes along each of the boundaries of the unit cell were 
coupled together to require the sides to remain straight and parallel to their original 
positions. This prevents the material being modeled from forming gaps along these sides. 
By applying a homogeneous normal stress in the longitudinal direction (o,,), the overall 
macroscopic, or effective, longitudinal elastic modulus (£„) and Poisson’s ratio (v , 2 ) can be 
determined by: 


E 


11 ■ 


( 1 ) 
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The transverse modulus (£,,) and minor Poisson's ratios (v,, and v^) can likewise be 
determined by applying a homogeneous normal stress in the transverse direction (a,,): 



V 


21 ' 



V 


23 * 



The requisite strain components are calculated from the appropriate 
components of the unit cell boundaries. 



(3) 


(4) 


C^) 

displacement 


Figure 2. Quarter-fiber unit cell for Haversian compact bone. 





Mori-Tanaka Method 

Zhao. Tandon and Weng (9j have recently applied the Mori-Tanaka method to the 
analysis of porous materials. The work represents an extension of previous developments 
for composite materials but the inclusions have been replaced by voids. The voids are 
considered to be spheroidal in shape, ranging from prolate to spherical to oblate. 
Expressions are presented for the effective constitutive constants of porous materials in 
terms of porosity, the aspect ratio of the voids, the properties of the solid phase, and the 
Eshelby transformation tensor. Several cases of assumed distributions and orientations have 
been analyzed. 

Perhaps the simplest case to consider as a model for Haversian compact bone is that 
of ‘needle-shaped’ voids, which is basically the same as a continuous fiber reinforced 
composite with voids instead of fibers. That is. the voids are long and slender, parallel to 
one another, and aligned with the overall longitudinal axis. The aspect ratio of the 
spheroidal voids is taken to be very high creating long extremely prolate voids. The 
macroscopic behavior is then transversely isotropic and the five independent constants are; 



(6a) 
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(6b) 
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The material properties of the isotropic solid phase are Eq, v q, Cq and ko- which represent 
the elastic modulus. Poisson’s ratio, shear modulus, and plane strain bulk modulus, 
respectively. The porosity is c, (volume fraction of voids) and the volume fraction of solid 
is Cfl. The effective properties are the longitudinal elastic modulus (En)> Poisson’s ratio 
(vjj), shear moduli (G,, and C,,). and plane strain bulk modulus (k,]). The transverse 
elastic modulus (Ej,) is readily calculated from these by; 
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Since osteons, and therefore Haversian canals, are not perfectly aligned with the longitudinal 
axis of the bones, it is also useful to transform the resulting effective moduli through a mild 
angle (0) to determine the longitudinal modulus when the axis of the voids is rotated off the 
main longitudinal axis. The transformed longitudinal modulus aligned with the global x-axis 
is then given by the familiar transformation equation; 
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cos*6 

sin^6 

^ 1 j 

JL- 




' 4 

.®12 ^11 , 


( 8 ) 


Another case examined as a model for compact bone is that of three-dimensional 
randomly oriented voids. The voids are again spheroidal but are considered to be randomly 
oriented in three-dimensional space creating a macroscopically isotropic material. This not 
only accounts for the non-alignment of Haversian canals but also for other sources of 
porosity such as Volkmanns canals and lacunae. Expressions for the effective shear 
modulus (G) and bulk modulus (k) are given in detail by Zhao. Tandon and VVeng [9], but 
the complexity prohibits repeating them here. The effective elastic modulus can then be 
calculated by: 


9kG 

3k*G 


(9) 


RESULTS AND DISCUSSION 

Predictions of the effective longitudinal modulus (£,,) using the finite element unit 
cell analysis are summarized in Fig. 3. The open squares are model results and the filled 
squares are experimental data from Schaffler and Burr (10). The heavy solid line is a 
power-law expression fit to the data by a least squares regression in log-log coon'.inates [ 10]. 
The correlation coefficient is 0.84 and the equation is; 

E = 33.88(1(in GPa) (10) 


The most striking feature of these results is that the experimental data exhibit much greater 
sensitivity to variations in porosity than do the model predictions. The finite element results 
decrease essentially linearly with porosity as would be expected for the unit cell 
configuration employed. The experimental results, however, decrease much more rapidly 
with increasing porosity, especially for the lower porosity specimens. 

A similar trend is also apparent in Fig. 4 for predictions from the Mori-Tanaka 
approach for needle-shaped voids. The modulus results are all normalized in this case, 
using a value of 33.88GPa from equation (10) for The solid curve labeled '0 deg’ is for 
the long cylindrical voids aligned with the global x-axis, in which case is given simply 
by equation (6a). This linear variation with porosity is basically the same as the finite 
element unit cell results. Rotating the axis of the voids off the global x-axis increases the 
sensitivity of the longitudinal modulus to variations in porosity but not dramatically. The 
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curves in Fig. 4 labeled 'IS deg’ and '30 deg' are actually £„ as given by equation (8) for 
0 of 15° and 30°, respectively. The longitudinal modulus decreases more because the 
transformation to the global j:-a.xis incorporates the effect of the modulus transverse to the 
axis of the voids, i.e. in equation (7). £„ is shown for reference in Fig. 4 as the ’90 deg’ 
curve. TTiis also serves to highlight the fact that even the most extreme case for the needle- 
shaped model (90°) does not produce as severe a variation in longitudinal modulus with 
porosity as reported experimentally. 

Putting too much emphasis on the results of Schaffler and Burr [ 10] as a basis for 
judging the usefulness of the analysis may. however, constitute an unfair comparison in 
several respects. First. Schaffler and Burr point out themselves that since their test 
specimens came from bovine bone tsteers). the lower porosity specimens likely represent 
mostly plexiform bone whereas the higher porosity specimens should contain more 
secondary osteonal (i.e. Haversian) bone. Plexiform bone has a microstructure quite 
different from Haversian none. The lamellar bone in plexiform bone is arranged in layers 
circumferentially around the medullary canal in the center of long bones in somewhat of a 
brick and mortar configuration (with elongated brieks). Porosity in plexiform bone comes 
primarily from vascular channels but they traverse in fairly irregular directions. Thus, it is 
likely that a high percentage of these vascular channels deviate considerably from the main 
longitudinal axis of the bone. .Another important consideration is that the ash content of 
the test specimens ranged from 0.664 to 0.723. Ash content is a measure of the degree of 
mineralization of the tissue. No statistically significant correlation was found between ash 
content and modulus for the 20 specimens, tjut Currey [11] has documented a strong positive 
correlation between modulus and calcium content for compact bone from a wide variety of 
species. The variation in ash content therefore almost certainly contributes to the variation 
in modulus for the experimental data plotted in Figs. 3 and 4. but the models do not account 
for this variable in their present form. It should be pointed out that the results of Currey 
(11) have not been used for comparison with the analysis herein because the wide variety 
of species makes it highly unlikely that much Haversian bone is included. Haversian bone 
is most prominent in adult humans but is mixed with plexiform and other bone types in 
many vertebrate animals. 

The Mori-Tanaka method for three-dimensionallv oriented voids would be expected 
to capture the more general spatial orientation of the vascular channels of plexiform bone 
as well as the inherent misalignment of Haversian canals in Haversian bone. The results 
shown in Fig. 5 seem to confirm this. The data and curve fit are the same as before, but 
the other curves represent the results of Zhao. Tandon and Weng [9] for several aspect 
ratios (a). The aspect ratio of 10 represents a prolate spheroidal void and the other three 
aspect ratios (.10, .05. .02) represent various oblate spheroidal voids. Prolate voids are more 
‘rod-like’ while the oblate voids are more disk-like'. The lower the aspect ratio for the 
oblate voids the flatter, or thinner, the disk. Fur oblate voids in general, and particularly 
so for very ‘thin disks’, proportionally more of the total volume of the void space is taken 
out of the plane perpendicular to the axis of loading. That is, oblate voids reduce the ‘net 
section’ more than do prolate voids for the same level of porosity. This gives rise to the 
more dramatic decrease in modulus with porosity for oblate voids. The dashed line curves 
in Fig. 5 for aspect ratios of .05 and .02 show the greatest sensitivity to variation with 
porosity and bound the experimental data. 

Improved predictions for strictly Haversian bone can be achieved through a 
somewhat ad hoc modification to the Mori-Tanaka method for needle-shaped voids. The 
properties given by equations (6a)-(6e) are retained to represent the basic configuration of 
long cylindrical voids aligned in the longitudinal direction. Equation (7) for £,„ however 
is replaced by an alternate expression for the modulus for uni-directionally aligned oblate 
spheroidal voids. The expression chosen is actually for the longitudinal modulus, or £„, in 
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Figure 5. Comparison of results from .\1ori-Tanaka method for 3-D randomly 
oriented voids with e.xperimental data (10). 


this case and is given in [9], The aligned oblate spheroids represent the effect of the 
Haversian canals, which run perpendicular to ihe a.\is of loading for the transverse modulus 
(£,,). Tlius. the modulus £,, for the uni-directionally aligned oblate spheroidal solution is 
used for £_._. in equation tb). Results are plotted in Fig. o for two aspect ratios i.02 & .01). 
The U deg' curves arc included for reierence. The oiner two curves from the analysis are 
for off-a.\is angles of 20° and 30°. Tlie oniv experimental data shown is for porosities of 
or greater. The assumption is that the results for higher porosities more likclv represent 
exclusively Haversian bone and not plexiform bone. Linear and power-law regressions of 
the 7 data points yield the following expressions; 

£ = 26.23-150.98£ (in G£a; r-= 80) 


£ = 28.18(1-£)*" {m GPa-, r-^Si) (Hb) 


which are plotted in Fig. 6 as well. The slopes of the '20 deg' and ’30 deg' curves for a = .02 
(Fig. 6(a.)) are at least similar to those of the curves fit to the data. Further decreasing the 
aspect ratio to .01 (Fig. 6(b.)) provides the best match between predictions and experiment, 
especially for an off-axis angle of 30°. It should again be emphasized, however, that no data 
currently exists in the literature on mechanical properties of exclusively Haversian compact 
bone as a function of porosity. While the present results are encouraging and increase the 
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prospects that theoretical models can eventually be devised to predict the mechanical 
properties of compact bone, more definitive development must be aecompanied by 
additional expertmental data. 


CLOSURE 

This comparison of micromechanics methods for modeling Haversian compact bone 
suggests that the Mori-Tanaka approach holds reasonable promise for predicting the 
variation in mechanical properties with porositv. The experimental data is scant, however, 
so the findings are quite tentative. The 2-D finite element unit cell analysis is limited in its 
present form to predicting an essentially linear decrease in modulus with porosity. The 
finite element method nevertheless posse.sses significant advantages: 11) several solid phases 
(including the cement line interphase material) can be readily incorporated. (2) the 
anisotropy of constituents can also be included easily, and (3) imperfect bonding (or even 
slip) between phases can be modeled as weil. Several options are therefore being 
investigated for improving the unit cell approach. These include extending it to three 
dimensions and/or devising alternate 2-D unit cell configurations. Current efforts are also 
considering possible ways to further improve and adapt the Mori-Tanaka method for this 
problem. The solid phase, or matrix material, in the current approach has isotropic material 
properties only, yet osteons and lamellar bone are anisotropic, at least to a certain extent. 
.■\nother approach being examined as this line of research continues is to combine the unit 
cell and Mori-Tanaka methods in some fashion to make use of the relative strengths of 
each. 
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ABSTRACT 

A finite element micromechenici model of Haversian compact 
bone has been used to investigate the influence of the cement line 
on the anisotropy of the macroscopic mechanical propenies. The 
micromechanics approach utilizes a 2-dimensional unit cell of 
material in a plane tiansverie to the main axis of the bone. The 
microstructure is idealized as consisting of three components: 
secondary osteons, cement line, and interstitial bone. Osteons 
are assumed to be arranged in a regular hexagonal amy parallel 
to one another. All constituents are assumed to be isotropic. 
The macroscopic mechanical properties examined are the 
longitudinal elastic modulus and the transverse elastic modulus. 
The ratio of the two is taken as an indicator of the degree of 
anisotropy. The thickness of the cement line has been varied 
from Ipm to fpm and its elastic modulus varied from 30% of 
that of the osteon to 1%. The effects on the longitudinal 
modulus are minimal, but the traiuvene modulus changes more 
significantly. For most of the range of parameters varied the 
modulus ratio is lower than experimentally observed, but the 
ratio approaches the experimental range at the extremes of higher 
cement line thickness and lower cement line modulus. 


INTRODUCTION 

Micromechanics modeling and analysis provides explicit and 
deterministic insight into the relationships between the 
mierostructure and macroscopic mechanical behavior of a variety 
of engineering materials. This approach typically entails 
i iealizing the geometry and architecture of the heterogeneous 
constituents of the material and identifying a repeating unit ceil 
of material that can then be analyzed in detail. A valid 
micromechanics model for compact bone tissue could be a 
powerful tool in helping elucidate the specific mechanical roles 
of the various microstnictunl elements. This would not only 
contribute directly to the basic science of bone tissue mechanics 


but it also r e p res ents an essential tool in the newly emerging field 
of biomimetics. Biomimetics refers to efforts to examine 
naturally occurring composite materials (e.g. wood, bone, insect 
exoskeletons, etc.) and identify salient 'design features’ that can 
then be mimicked in producing improved man-madc.composiic 
materials. 

Initial efforts to model bone as a composite material treated it 
as a two-phase material of collagen and apatite crystals (Currey, 
1964, 1969; Pickarski, 1973). Katz (1976, 1981) proposed a 
hierarchical model for compact bone with the 'microstructural' 
level being that for which secondary osteons, or Haversian 
systems, could be considered to be the main microstructural 
element. On this level, the osteons are considered to be 
analogous to the fibers of a fiber-reinforced composite material 
(although hollow) and the cement line and interstitial bone are 
taken together as a binder material forming the matrix of the 
composite. The same basic approach was followed in developing 
the current micromechanics model (Hogan, 1988, 1990, 1992). 
which utilizes a finite element analysis of the material unit cell to 
permit more general and explicit treatment of microstructural 
details. Others have also used a similar approach in studying 
Haversian compact bone (Crolet, 1990; Pidapartidb Burr. 1991). 

The objective of the present study has been to use the 
micromechanics model to examine the potential role of the 
cement line in determining the anisotropy of the mechanical 
properties of Haversian compact bone. A prominent ai^ 
consistent cement line surrounding secondary osteons >* 
distinctive of Haversian compact bone, and its meehanical role 
has been the subject of much study and sp ecul ati on . The 
macroscopic mechanical properties of compact bone ate tstba 
well documented to be anisotropic, either transversely isotropic 
or ortbotropic. In either case, the longitudinal cl a s tic modulus is 
typically much greater than the elastic modulus transverse to the 
longitudinal axis. The potential sources of this anisotropy mclu^ 
porosity and anisotropy of the osteons or interstitial bone, but Iha 
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FIGURE 1. IDEALIZED MICROSTRUCTURE OF 
HAVERSIAN COMPACT BONE. 


purpose of (he present study is to eumtne (he extent to which the 
anisotropy can be attributed to the cement line and its 
characteristics. 

MICROMECHANICS MOOEUNG 
The idealized microstructure of Haversian compact bone tissue 
is shown in Fig. 1. The material is assumed to extend 
continuously perpendicular to the transverse plane shown 
allowing a Z-dimensional micromechanics model. The only void 
space included explicitly in this model is the Haversian canal. 
Three components make up (he solid phase; (1) secondary 
osteons representing a hollow Hber component, (2) interstitial 
bone representing the matrix material, and (3) the cement line 
representing a distinct fiber/matrix ’interphase' material. The 
cement line is assumed to be perfectly bonded to the osteon and 
interstitial bone. The osteon fibers are spaced in a regular 
hexagonal array. The rectangular unit cell has been identified 
based upon symmetry of the material microstructure as well as 
symmetry of loading conditions. 

A typical finite element mesh for the model is shown in Fig. 
2. The model consisted of 322 8-noded isoparametric elements 
and 10S3 nodes. The elements include a generalized plane strain 
formulation to allow out-of-plane loading perpendincular to the 
plane of the model. All nodes on the straight boundaries of the 
unit cell are constrained to deflect such that the boundaries 
remain straight and parallel to their original positions. This 
prevents the material from forming new internal 'voids' during 
deformation. The overall dimensions of the unit cell are lOOnm 
173.2nm. The outer radius of the osteon is ST.SSnm creating 



a fiber volume percent (or percent Haversian area) of 70% for all 
cases studied. The radius of the Haversian canal for the mesh 
shown IS 18.2iim yielding a porosity of 3%. A porosity of 5% 
was also studied by increasuig this radius to 23.5iim. The 
thickness of the cement line was either Ifim, 3Mm. or 5nm. All 
three solid components were modeled as linearly elastic and 
isotropic. The clastic modulus was taken to be 12GPa for the 
osteon and i5GPa for the interstitial bone. A Poisson's ratio of 
0.3 was used for both of these as well. The cement line was 
considered to be much more compliant than the bone 
components. Its precise mechanical properties are unknown, 
however, so the following values were used for the elastic 
modulus of the cement line: 6. 3.6, 1.2, 0.6, and 0.12 (all in 
GPa). As a percent of the osteon modulus, these correspond to 
50%. 30%. 10%, 5%, and 1%, respectively. A Poisson’s ratio 
of 0.4 was used for the cemem line. The model was loaded with 
a uniform applied stress in the horizontal direction to determine 
the transverse elastic modulus. The in-piane Poisson’s ratio can 
also be determined fiem this loading case using the lateral 
contraction displacements of the upper and lower boundaries of 
the unit cell. The model can also be loaded in the longitudinal 
direction (perpendicular to the plane of the model) to determine 
(he longitudinal elastic modulus and additional Poisson's ratios. 
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TABLE 1 

MACROSCOPIC LONGITUDINAL MODULUS 


TABLE 3 

LONGITUDINAL TO TRANSVERSE MODULUS RATIO 


Ed. 

(%l 

1//m 

3% 5% 

3^/01 

3% 5% 

5/jm 

3% 5% 

50 

12.40 

12.16 

12.10 

11.86 

11.80 

11.56 

30 

12.36 

12.12 

11.99 

11.75 

11.61 

11.37 

10 

12.32 

12.08 

11.87 

11.63 

11.40 

11.17 

5 

12.31 

12.07 

11.84 

11.60 

11.36 

11.12 

1 

12.30 

12.06 

11.82 

11.58 

11.32 

11.08 


units in GPa 


TABLE 2 

MACROSCOPIC TRANSVERSE MODULUS 


Ed. 

(%) 

Tpm 

3% 5% 

3pm 

3% 5% 

5pm 

3% 5% 

so 

11.60 

10.98 

11.26 

10.65 

10.93 

10.33 

30 

11.46 

10.84 

10.87 

10.27 

10.33 

9.75 

10 

10.87 

10.27 

9.54 

8.99 

8.53 

8.05 

5 

10.18 

9.60 

8.31 

7.84 

7.09 

6.71 

1 

7.43 

7.04 

4.88 

4.68 

3.68 

3.55 


units in GPa 


RESULTS AND DISCUSSION 
The macroscopic elastic modulus predicted by the model is 
essentially the same as that given by a rule-of-mixtures (uniform 
strain) approach for the longitudinal direction. The cement line 
is a ver)' small percentage of the total solid phase componenu 
(ranging from 1.7% to 8.6%), so the variation in the longitudinal 
modulus with cement line properties and thickness is minimal as 
shown in Table 1. The cement line modulus is given as a 
percent of the osteon modulus. Results are presented for cement 
line thicknesses of l^m, 3^m, andSpra, and for porosities of 3% 
and 5%. As would be expected, the modulus decreases for lower 
cement line moduli and for higher cement line thicknesses. Note 
as well that the moduli are generally slightly lower for 5% 
porosity versus 3%. Aaisotropy in the macroscopic moduli is 
clearly evident in comparing the macroscopic transverse moduli 
in Table 2 with the corresponding longitudinal moduli in Table 
1. The transverse moduli are all lower than the corresponding 
longitudinal moduli. The results in Table 2 also indicate that the 
traiuverse moduli ate more sensitive to changes in cement line 
modulus and thickness. The transverse modulus drops off much 
more significantly as the cement line modulus decreases to 10% 
or less of the osteon modulus. At the extreme case of 1 % the 
tratuverse modulus ranges from 64% to 34% of the 
corresponding value for a cement line modulus of S0% of the 


Ed. 

(%) 

U 

3% 

rm 

5% 

3pm 

3% 5% 

3% 5% 

SO 

1.07 

1.11 

1.07 

1.11 

1.08 

1.12 

30 

1.10 

1.12 

1.10 

1.14 

1.13 

1.17 

10 

1.13 

1.18 

1.24 

1.29 

1.34 

1.39 

5 

1.21- 

1.26 

1.42 

1.48 

1.60 

1.66 

1 

1.66 

1.71 

2.42 

2.47 

3.09 

3.12 



Cement Line Modulus (% of osteon) 

FIGURE 3. MODULUS RATIO AS A FUNCTION 
OF CEMENT LINE PROPERTIES. 


osteon modulus. The effect is again slightly greater for the 5% 
porosity results. 

As a q u a n ti t ative indication of the degree of anisotropy in the 
macroscopic modulus predictions, the ratio of the longitudinal 
modulus to the tratuverse modulus has also been c a l c ul a ted . 
These resulu are presented in detail in Table 3 and summarized 
in Fig. 3 for S% porosity and the two extreme eemem line 
thicknesses only (l>im and Sum). The longitudinal modulus is 
only 10 to 20 percent higher than the traruverse modulus for 
cement line moduli of S0% and 30% of the osteon modulus. As 
the cement line modulus further decreases, however, the ratio 
increases rather dramatically, particularly for the thickest (Sum) 
cement line. Overall, these ratios are lower than those typically 
reported in the literature. Cowin (1989) summarized resulu for 
human compact bone from several sources, and these yield the 
following ratios: 1.48, 1.46, 2.66. 2.16, 1.67, and 1.49. 

At least two views are possible when comparing the modeling 
resulu with experimental dau. One is that the model simply 
lacks sufficient realism in iu presem form because k predku a 
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much lower moouius ratio for most of the range of parameters 
studied. In this case, other sources of anisotropy may need to be 
included in the model. Prominent among these would be to 
consider the osteon to be anisotropic due to the preferential 
onenuuon of collagen and apatite in the concentric layers of 
osteonic lamellar bone. This can be modeled explicitly within 
the framework of the present approach but the precise degree of 
anisotropy of the osteon is unknown. This becomes to a large 
extent then another adjustable parameter that can be varied over 
a wide range as needed. Osteon anisotropy has been mvestigated 
theoretically (Katz. 1981; Pidaparu et ah, 1992), but no direct 
experimental evidence is currently available to firmly establish 
anisotropic osteon properties. Another view is that the lower 
cement line modulus cases in the present analysis more accurately 
reflect the true nature of the cement line and this is why the 
predictions match the experimental data more closely for these 
cases. The precise mechanical properties of the cement line are 
also unknown at this pomt, but most consider it to be very 
compliant (Burr et ah, 1988; Frasca, 1981; Schaffler et ah. 
1987) and perhaps even viscous (Lakes and Saha. 1979). The 
current model treats the cement line as perfectly bonded to the 
surrounding bone, which may therefore overstiffen the model 
relative to the actual physical situation. The lowest cement 'ine 
modulus values of S% and I7c could be considered somewhat 
unrealistic since they are 20 and 100 times, respectively, lower 
than the osteon modulus, but the net effect taken together with 
the perfect bonding may in fact combine to produce a response 
quite close to the actual case, in either case, the overriding need 
is for more definitive experimental data on the mechanical 
propeities of the osteon and the cement line. Such information 
will not be easy to come by, however, and the present study 
demonstrates the uttlity of micromechanics modeling in exploring 
the possible effects of microconstituent properties on macroscopic 
properties. Specifically, the results herein suggest that the 
cement line is likely a significant contributor to Haversian 
compact bone anisotropy. The macroscopic anisotropy in 
Haversian bone mechanical properties is not necessarily due 
solely to osteon anisotropy and can conceivably be explained 
based upon cement line properties alone. These possibilities 
should therefore remain open to mquiry until direct experimental 
evidence mdicates otherwise. 
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THEORETICAL MODELING OF THE EFFECTS OF POROSITY ON 
THE MECHANICAL PROPERTIES OF COMPACT BONE TISSUE 

H. A. Hogan and R. J. DeFrese 

Dept, of Mechanical Engr., Texas A&M University, College Station, Texas, USA 
INTRODUCTION 

At the microstructural level compact bone contains a network of canals and cavities that house 
blood vessels, nerves, and cells. The porosity of these void spaces can range from 2% to 20%. A 
basic first step in establishing theoretical structure/property relationships for compact bone tissue is 
to account for the separate effects of void space and the solid phase. Compact bone tissue can 
actually be classified into several different types depending on the details of the microstructure. The 
specific type of interest in the present work is Haversian compact bone tissue, which has secondary 
osteons as the predominant microstructural element. These osteons are commonly considered to be 
approximately cylindrical in shape and aligned roughly with the long axis of the bone. The main 
void space in Haversian compact bone is the Haversian canal, which is a hollow cylindrical canal 
that runs down the center of each osteon. The objective of the present study has been to evaluate 
and assess various theoretical approaches for modeling these types of void spaces and their effects 
on mechanical properties. As a beginning point, the mechanical property of interest is simply the 
elastic modulus of the material. 

POROSITY MODELING 

Most attempts to develop analytical expressions relating porosity and mechanical properties 
of various engineering materials have resulted in empirical correlations based upon extensive 
experimental data. These expressions typically take the form of power law relations with particular 
coefficients and exponents for particular materials. Theoretical (Zhao et al., 1989) and semi- 
empirical (Bert, 1985) treatments are also available and have been examined as candidates for 
modeling the cylindrical voids of Haversian compact bone tissue. Zhao et al. (1989) use an effective 
medium approach to develop expressions for the elastic constants of materials with spheroidal voids. 
The voids can be either aligned or randomly oriented and can have aspect ratios less than one 
(oblate) or greater than one (prolate). Cylindrical voids are a special case in which the voids are 
aligned and have an aspect ratio approaching infinity. Bert (1985) developed a combined theoretical 
and empirical treatment for cylindrical and spheroidal voids. The cylindrical and spheroidal voids 
can be taken to be arranged in a square array or a hexagonal array. Finite element micromechanics 
modeling has also been used for studying compact bone tissue (Hogan, 1992). Porosity due to the 
Haversian canal is modeled explicitly, and distinct osteon, cement line, and interstitial bone 
components are included separately for the solid phase. 

RESULTS AND DISCUSSION 

Two fundamental factors to be considered in a theoretical model for porosity effects in 
compact bone are; (1) the variation in mechanical properties with porosity; and (2) the degree of 
anisotropy in the mechanical properties due to the presence of porosity. Elastic modulus predictions 
are compared in Fig. 1 for the analytical models of Bert (1985) and Zhao et al. (1989) and for the 
current finite element micromechanics model. The moduli values are normalized by dividing by the 
elastic modulus of the sol'd phase. All three approaches yield the same results for the elastic 
modulus in the axial direction, i.e. along the main longitudinal axis of the long bone. The cylindrical 
voids are also aligned with this axis so the decrease in modulus is linear and proportional to the 
porosity. The elastic modulus transverse to the longitudinal axis is more sensitive to changes in 





porosity. The results for Bert's model are for a hexagonal array of cylindrical voids and show a 
slightly greater decrease in modulus with porosity than the other two approaches. The rate of 
decrease, however, is greatest for the finite element results. The degree of anisotropy predicted by 
the models is depicted by plotting the ratio of the axial modulus to the transverse modulus. Bert's 
model predicts greater anisotropy due to the lower transverse modulus values. The ratio increases 
from approximately 1.0 to 1.35 for Zhao et al. and the finite element model, but reaches a maximum 
of 1.55 for Bert. These are generally lower than typical ratios reported for human compact bone. 
Cowin (1989) summarized data from four sources which yield the following ratios: 1.48, 1.46, 
2.66, 2.16, 1.67, 1.49. Six values are included because two of the studies found the bone to be 
orthotropic so two different transverse moduli were reported. None of the studies included any 
microstructural information about the specimens tested so the porosities are not known. An 
important point to consider as well is that the two analytical approaches assume the solid phase to 
be homogeneous and isotropic. Each component of the solid phase in the finite element model is 
likewise treated as isotropic. The observed degree of anisotropy is likely due in part to the 
microstructural architecture of the solid phase as well as the presence of voids. The theoretical 
models nevertheless clearly demonstrate that a significant portion of the anisotropy of Haversian 
compact bone can be attributable to porosity. 

Only limited experimental data is available in i s 
the literature on the variation in elastic 
modulus of compact bone with porosity. 

Currey (1988) examined the influence of 
porosity and mineralization on the elastic ^ 

modulus of bone 18 different species of J 12 

animals but none were from human sources. ^ „ 

Volume fraction (1-p) was correlated using a 
power law expression with the exponent J o.s 
ranging from 3.1 to 3.5. Schaffler and Burr g 
(1988) report a power law exponent of 10.92 o ° ® 
for modulus as a function of volume fraction 0.4 
10.92. These are not linear relations as the 
models herein predict, but neither study 
included human bone, which is predominantly 
Haver-sian bone. Other species are typically . 

a mix of plexiform and Haversian bone, so the >''8 ' '’"‘‘•ebons of modulus vs porosity. 

models would not be expected to predict well for these. 
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